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The large-scale/long-time transport of inertial particles of arbitrary mass density under 
gravity is investigated by means of a formal multiple-scale perturbative expansion in 
the scale-separation parametre between the carrier flow and the particle concentration 
field. The resulting large-scale equation for the particle concentration is determined, 
and is found to be diffusive with a positive-definite eddy diffusivity. The calculation 
of the latter tensor is reduced to the resolution of an auxiliary differential problem, 
consisting of a coupled set of two differential equations in a (6-|-l)-dimensional coordinate 
system (3 space coordinates plus 3 velocity coordinates plus time). Although expensive, 
numerical methods can be exploited to obtain the eddy diffusivity, for any desirable 
non-perturbative limit (e.g. arbitrary Stokes and Froude numbers). The aforementioned 
large-scale equation is then specialized to deal with two different relevant perturbative 
limits: i) vanishing of both Stokes time and sedimenting particle velocity; ii) vanishing 
Stokes time and finite sedimenting particle velocity. Both asymptotics lead to a greatly 
simplified auxiliary differential problem, now involving only space coordinates and thus 
easy to be tackled by standard numerical techniques. Explicit, exact expressions for the 
eddy diffusivities have been calculated, for both asymptotics, for the class of parallel 
flows, both static and time-dependent. This allows us to investigate analytically the role 
of gravity and inertia on the diffusion process by varying relevant features of the carrier 
flow, as e.g. the form of its temporal correlation function. Our results exclude a universal 
role played by gravity and inertia on the diffusive behaviour: regimes of both enhanced 
and reduced diffusion may exist, depending on the detailed structure of the carrier flow. 



1. Introduction 

Transport of passive particles is a problem of major importance in a variety of do- 
mains r anging from astrophysics a nd geophysics to technological applications and biol- 
ogy (.Ruiz. Macias. fc Petersll2004 ). A quantity of particular interest is the rate at which 
particles are transported by the flow. 

For inertialess particles (i.e. fluid particles), for times large compared to those char- 
acteristic of the given velocity field, trans port is usually d iffusive and characterized by 
effective (enhanced) diffusion coefficients (jKraichnan 19871) . the so-called eddy diffusivi- 
ties. These coefficients incorporate all the dynamical effects played by the velocity field 
on the particle transport. 
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Although the diffusive scenario with eddy diffusivities is the typical one, there exist 
cases where superdiffus ion is observed alrea dy for simple incompr essible laminar flows 
( Castiglione et al. 1999f) and synthetic flows ( Andersen et al. 20001 ). 
From an applicative point of view, eddy diffusivities have long been fruitful concepts 
in turbulent-transport theory, and their use has made the computation of turbulent- 
transport problems_20Ssible at Peclet numbers too high for full numerical simulations 



(see, e.g.. lLesieuij|l997l ). 



In an asymptotic case of interaction between modes whose space and time scales are 
strongly separated, eddy diffusiviti es have been ca l culate d exploiting a multiple-scale ex- 
pansion in the scale parametre bv iBiferale et al. (1995). Conversely, for realistic flows, 
active on all space-time scales, eddy diffusivities are gener ally dependent o n all flow 
characteristics and no general expression for them is known ( Kraichnanl 1987 ) . 

The sit uation is even more complex i n the case of transport of particles wit h inertia 

(lGatignoll983; Max ey fc Rjlcy 1983- Balkovskv. Falkovich fc Fouxon 2001; Wilkinson fc Mehhe 
2003t ICencini et aLH2006t IVolk et aLlboOSf ). The main difficulty arises from the fact that 



in this case the fluid velocity does not coincide with the particle velocity, with the result 
that a larger phase space has to be considered (i.e. involving both particle position and 
velocity). The standard Fokker-Planck equation (which only involves space variables), 
holding for the concentration fleld in the case of fluid particles, is replaced, for particle 
with inertia, by a Kramers equation for the phase-space density now involving both space 
and velocity coordinates. 

In many cases of interest, the dependence of the phase-space particle density on the 
velocity coordinates is not of particular interest (a relevant example is in the realm of 
pollutants dispersion, where the crucial objective is to predict the space-time behaviour 
of the pollutant concentration rather than the velocity distribution of the pollutant) and 
can be averaged out in a way such that only space-time variables are involved. 
Closed eq uations for such an obser vable (the so-called marginal density) have been ob- 
tained bv IPavliotis fc Stuart ( 2005 ) exploiting a formal multiple-scale expansion in the 
scale-separation parametre. The attention was focused there on the sole case of heavy 
particles and in the absence of gravity. The latter is however relevant in many situations 
of interest, including those where the particle acceleration is (at least) of the same or- 



der of magnitude of the gravitational acceleration (see, e.g.. iMaxevI 119871 ). This latter 
case corresponds to flnitc terminal sedimenting velocity even in the presence of particles 
with small inertia. A question which naturally arises is on whether a constant (gravity- 
induced) falling/rising terminal velocity is accompanied by a correction to the diffusion 
process occurring in the absence of gravity. Answering this question is one of the main 
contribution of the present paper. 

To achieve our aim we have exploited a formal multiple-scale expansion in the scale- 
separation parametre on the Kramers equation for the particle probability density func- 
tion. Different particle densities have been considered, ranging from the limit of heavy 
particles to the one of light particles (i.e. b ubbles) . The spirit of the work is similar to the 
one which motivated Biferale et al. ([l99i). Unfortunately, for the reasons already men- 
tioned, the analytical treatment is here much more cumbersome, a fact which imposes 
a more formal treatment of the multiple-scale expansion. With the hope of making the 
present manuscript more readable, we decided to postpone the most technical issues in 
dedicated appendices. Here is a short summary of our main achievements. 
We found that the large-scale/long-time behaviour of the particle concentration field is 
diffusive and the eddy diffusivity is positive defined. This holds for all Stokes times and 
particles density. To obtain the eddy-diffusivity tensor one has to solve (numerically) a 
system of two, one-way coupled, partial differential equations (the auxiliary problem) in 
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a phase space involving both the particle coordinates and their velocity (plus time). In 
a three-dimensional space, the differential problem thus involves 7 coordinates. Conse- 
quently, the numerical solution requires a considerable computational effort. 
In order to simplify the problem, and at the same time motivated by situations of inter- 
ests where the Stokes time is small, a further perturbative analysis has been carried out 
in the limit of small inertia. Two asymptotics have been considered: one corresponding 
to small Stokes and vanishing terminal sedimenting veloci ty, the o ther corresponding to 
small Stokes and finite terminal sedimenting velocity as in lMaxevl (,1987). In both cases, 
the resulting expressions for the eddy-diffusivity tensor greatly simplify: they still involve 
an auxiliary problem, but now in the three-dimensional space. Also in this case, in gen- 
eral, only numerical approaches are capable of obtaining the eddy-diffusivity tensor. At 
the perturbative order here considered, gravity does not affect the diffusion process in 
the first asymptotics, while it does in the second one. 

As we ment ioned above, the eddy diffusivity is a very old idea in turbulence (see, e.g., 
iFrisch 1995f) which permitted to simulate high Reynolds/Peclet number flows. Motivated 
by this consideration, we have isolated from both our asymptotics (vanishing or finite 
sedimenting velocity) a class of relevant flow fields from which the eddy diffusivity can be 
computed analytically. This is indeed the case for the class of parallel flows. The aim was 
to understand if gravity plays to enhance or reduce the rate of transport. Answering this 
question might help researchers involved in closure problems to try to parameterize the 
role of gravity in the diffusion process in terms of relevant properties of the flow flcld. Our 
results are not encouraging from this point of view: we performed a detailed sensitivity 
study, from which it clearly emerges a dependence of the effect played by gravity and 
inertia on the details of the carrying flow. A simple, universal, way to parameterize the 
role of gravity and inertia on the diffusion process seems to be very difficult to achieve. 

The paper is organized as follows. In § [2] we sketch the problem under consideration 
by recalling the significant equations, and we perform a multiscale expansion to obtain 
the eddy-diffusivity tensor and the auxiliary equation. In § [3] we analyse our small- 
inertia expansion order by order. Section |4] shows explicit results for parallel flows and 
namely for a speciflc case, the well-known Kolmogorov flow (steady or random in time). 
Conclusions and perspectives follow in § [Sj The appendices § [X] and § |B] are devoted to 
show the mathematical details of the calculation and to recall some basic notions about 
the Ornstein-Uhlenbeck process, respectively. 



2. Large-scale dynamics and multiscale expansion 

The model of transport we consider here refers to point-like inertial particles subject 
to a constant (gravitational) acceleration g (for the sake of generality, in d spatial dimen- 
sions) and carried by an incompressible velocity fleld u{x, t). Beside incompressibility, we 
suppose that the velocity field is steady or periodic in time (with period T) , and periodic 
in space with unit cell B of linear size t. It is not a restriction to focus on velocity fields 
the average of which vanishes over B: 



/ da; u{x^ t) — 



(2.1) 



The same technique can be extended to handle the case of a random, homogeneous, 
and stationa ry velocity field with some n ontrivial modifications in the rigorous proofs of 
convergence (lAvellaneda fc Maidalll99lh . 

Neglecting any feedback effect on the transpor ting fluid and taking only into account 
the added-mass effect in a simplified way as in iMartins Afonsol (|2008f ). the Lagrangian 
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dynamics reduces to the foUowing set of stochastic differential e quations for particle 
position iX{t)) and co-vclocity (V(t)) (|Maxev fc Rilevlll983l : lGatignoI,1983i) : 



X{t) = V{t)+l3u{X{t),t) 

V{t)~ {l-l3)u[Xit),t] 



+ (1 - P)9 + 



/2k 



(2.2) 



In (|2.2p ri{t) is a white- noise process coupled by a constant difFusivity k (|Reekslll988h . 
the pure number /3 is the added-mass factor f3 := 3pf/(pf + 2pp) £ [0,3] built from 
the constant fluid (pf) and particle (pp) mass densities. Non-vanishing values of /3, i.e. 
particles not much heavier than th e fluid, in duce a discrepancy between the particle 
velocity, X{t), and co- velocity, V{t) ( Bedl2003l ) . Finally, the Stokes time rs expresses the 
typical response delay of particles to flow variations. For spherical inertial particles of 
radius R, the Stokes time is related to the kinematic viscosity v of the carrier fluid by 

rs - RyiSi^P). 

A s tandard application of stochastic calculus (see, e.g. . IChandrasekhan 1 943t I Van Kampenl 
2007t) yields the Kramers equation for p{x,v,t), the phase-space density: 



d_ 

dt 



d_ 
dv 



+ (l-/?)9 



p = 0. (2.3) 



Our aim is to derive the large-scale asymptotics of the marginal density P, 



P{x,t) i dv p{x,v,t) , 



(2.4) 



describing variations of the configuration-space marginal distribution over spatial and 
temporal scales 0{x) ~ L ^ £ and 0{t) ~ T :s> T, respectively. The technical tool 
we employ to extricate this large-scale asymp totics from the full dynamics is known 
as homoqenization or multiple-scale expansion ( | Bensoussan. Lions fc Papanicolaoulll978 : 
Bender fc Orszad 19781 Pavliotis fc Stuart 2007 ). The idea is to use the scale-separation 
parametre 

e := i/L < 1 

to introduce a new set of ( "slow" or "large-scale" ) variables 



(2.5) 



which we shall treat as independent from the corresponding "fast" variables x and t. 
We do not introduce any slow variable V associated to the co- velocity v, since we are 
interested only in the asymptotics of the marginal configuration-space distribution (12. 4p . 
We define instead the larg e-scale time variables T and to de couple diffusive from 
ballistic degrees of freedom (jMazzino. Musacchio fc Vulpianill2005l ). We also hypothesize 
that the expression of the initial conditions for p.3p in terms of fast and slow variables 
takes the form 

p{x, X, V, 0, 0, 0) = p^ix) p,{v) P{X) , (2.6) 

where px has the same periodicity of the velocity field u. The product form (|2.6p is 
non-restrictive owing to the linearity of the dynamics (|2.3p . This latter is couched by the 
change of variables p.Sp , 




into the form 



with 
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+ (1-/3)9 



(2.7) 



- — A„ , (2.8a) 



(2.8&) 



r(2) 



(2.8c) 



Our program is then to solve (j2.7l) in a power-series expansion in the scale-separation 
parametre, 



p{x, X, V, t, T, T*) = ^ (a;, X, v, t, T,T^) , 

with initial conditions 

(a;, X, V, 0, 0, 0) = SnopAx)Pv{v)P{X) , 
based on the "unperturbed" Cauchy problem in the fast variables, 

£(°)p(°) = , 

p^°\x, X, V, 0, T, Tt) = pAx)pAv)P("^ {X, T, T*) . 



(2.9a) 

(2.96) 

(2.10a) 
(2.106) 



The final objective is to attain closed equations for the large-scale dynamics, at the lead- 
ing order described by P'^'^\ in the form of solvability conditions for the perturbative 
hierarchy. As a result, the details of the small-scale dynamics are summarized in the def- 
initi on of renormalized transport coeffic ients, such as the effective terminal velocity (see 
also axevlll987l : iMartins Afonsoll2008l . and references therein) and the eddy diffusivity, 
for the advective and diffusive rescalings respectively. 

For what concerns the fast variables, here the aim is to find square-integrable solutions 
which enjoy the same space-time periodicity as the flow fleld u, and which decay suffi- 
ciently fast at infinity in the co- velocity varia ble. Under physically reasonable r egularity 
assumptions on u, the results of chapter 3 of 'BensijiKsari, Lions & Papanicolaou (1978) 



carry over to the operator L'^^\ Its adjoint on the space L^(I 



[o,r]), 



(1 - I3)u-v 



+ (1-/3)9 



|--4a., (2.11) 



has as zero modes only the constants with respect to the fast variables. By the Fredholm 

alternative, the ex istence and uniqueness of the inverse of ( Bensoussan. Lions &: Papanicolaou 
1978t Pavliotis fc S tuart 2007) is guaranteed by imposing that the kernel of £^°''t be or- 
thogonal to the perturbations spawned by (I2.86P and (|2.8cp . or in other words that the 
integrals of the right-hand sides of (|2.16p and ()2.23p be vanishing. This will allow us to 
determine the dependence on the slow variables. 
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2.1. Derivation of the large-scale evolution 
Since the Kramers evolution (12.31) is measure preserving, it is not restrictive to posit 

/ ^ [dx [ dv [ dX p{x,X,v,t,T,T^) = 1 . (2.12) 

Combined with (|2.9&p . ()2.12p allows us to fix the normalization of the terms of the 
multiscale expansion to 

[dx [ dv [ dX/,(") = 5„o . (2.13) 

Jo ' Jm JR'' Js.'' 

Since the solution of the unperturbed Cauchy problem (|2.10a| - [2.106p takes the form 

p^'\x,X,v,t,T,Ti)^p(^'\x,v,t)p("\X,T,Ti) , (2.14) 

with 

-0 , (2.15) 

from (|2.13p we obtain 

fdx f dvp'^"Hx,v,t)^l^ f dXpW(X,T,r*) . 
Jo ' Jb Jw^ JR<i 

We will show below that this condition is mathematically self-consistent because the 
solvability conditions yield for P a measure-preserving evolution law. Physically, this is 
desirable as it corresponds to the statement that, when focusing on the sole large scales of 
the problem, P can be interpreted itself as a mass density. It is also a necessary condition 
to correctly reproduce the well-known small-inertia limit for the eddy difFusivity of tracer 
particles. 

2.1.1. Solvability condition for p'^^h large-scale ballistic motion 
At O(e^) we have: 

d d 



('^ + ^"^ dX ' dTt 



. (2.16) 



By the Fredholm alternative the problem admits a unique solution if and only if 

[ ^ /da; / dt;/:Wp(o) = -/ ^ [ dx [ dv C^'>^ p^'^ = . (2.17) 
Jo I Jb Jr.'' Jo ' J« Jri^ 

Upon inserting (|2.14p into (|2.17p we get the solvability condition 

_|_p(0) ^x,T,Ti)+w Ap{0) (X, T, ) = , (2.18) 
with w the terminal velocity ( Maxevlll987 : Martins Afonsoll2008l ) specified by 

w -.^ [ ^[dxf dv [v + Pu{x,t)]p''°\x,v,t) . (2.19) 
Jo I Jb Jri^ 

The physical interpretation of p.lSp is to associate the terminal velocity to ballistic 
large-scale motion. The solution of (|2.18l) requires P^^-* to satisfy 

P^°\X,T,Ti) = P'^^^X -wT\T,0) 

wi thout fixing the functional d ependence with respect to the first two arguments. As 
in Mazzino. Musacchio &: Vulpiani (2005h we can leverage this fact to find the function 
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P*^\x,T) satisfying the Fredholm ahernative condition at the second order in e. We 
can then reconstruct the fuU solution by the Galilean transformation 

p(°)(x,T,r*) ■.^p!:°\x-wT^,t) . 

2.1.2. Solution of the equation for the leading-order correction p^^^ 
Exploiting (|2.18l) . the first-order equation (|2.16p reduces to: 

dp(o) 



C(o)pW =-p(") {v + /3u-w) 



dX 



(2.20) 



The initial condition (j2.9&P rules out the presence in p*-^^ of any solution of the associated 
homogeneous problem. We are therefore entitled to couch the solution into the form 



{x, X, V, t, T, Ti) = {x, V, t) ■ (X, T, T^) 

with d-dimensional vector field satisfying 

C^°^p'^^\x,v,t) = - [v + I3u{x,t) -w]p'^°\x,v,t) , 

p^'^\x,v,0) = . 



(2.21) 



(2.22a) 
(2.226) 



Equation ()2.22ap is usually called "auxiliary equation" or "cell problem" . The advantage 
of (|2.21l) is that it satisfies at glance the normalization condition (I2.13|) with n = 1, in 
consequence of the gradient form in the slow variable X . 

2.1.3. Solvability condition for p^^h large-scale diffusive motion 
Large-scale diffusive effects set in only at order O(e^). Namely, 

d d 



admits a unique solution if and only if 

dv (£(i)pW+£(2)p(0) 



{v + l3u) 



dX dTt 



(2.23) 



T 



T 



dt;£(°V^'^ -0. (2.24) 



Upon inserting on the right-hand side of (|2.24l) the explicit expressions (j2.8&P of ^C^^-* and 
(j2.21|) of p^^\ after straightforward algebra we obtain 



Ap(X,T,T^)=i.,.^£^P(X,T,T^), 



(2.25) 



where we have introduced the eddy-diffusivity tensor 



Ax 



Av {v -\- /3u{x,t) - w]^pI}^ {x,v,t) 



(2.26) 



with S the trace-preserving index symmetrization operator: SA 
However, finding explicit expressions for (|2.26p requires the knowledge of the terminal 
velocity w (|2.19|) and of the vector field p'^-* (|2.22al) . and thus also of the density p^^^ 
(|2.15p . These two latter partial differential equations are defined in a (2(i-f l)-dimensional 
space according to the definition of (|2.8ap . and are thus prohibitive to solve also 
numerically. In § |3] we will turn to the derivation of systematic approximations for these 
quantities in the limit of small inertia, which allows us to decouple the co- velocity degree 
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of freedom and therefore to recast the problem in terms of "innocent" forced advection- 
diffusion equations, which are {d+ l)-dimensional and can be solved relatively easily for 
specific flows (§ 2]) or in general from the numerical point of view. 

In order to prove the diffusive nature of the large-scale motion described by (|2.25p . it 
remains to show that the eddy-diff usivity tensor (12.261) is p ositive definite. The proof 



proceeds along the same lines as in iPavliotis fc Stuart (l2005h . therefore we only briefly 



recall its outline. By (|2.19l) and the Fredholm alternative, the auxiliary equation 

v,t)^v + Pu{x, t)-w (2.27) 

admits a unique solution. We can therefore exploit it to write (|2.26l) as a matrix element 
of An integration by parts and the use of (I2.22al) and of (|2.27p again, finally give 

K^.^sT^ [dx [ dvp^°^x,v,t)[v + puix,t) - w] C^"^^-' [v + I3u{x,t) - w]^ , 

which is essentially a weighted average of the projection of a positive-definite operator 
onto the complement of its kernel. 

The solution of (|2.25p satisfying the initial condition 

P^°\X,0) = P{X) 

yields the large-scale asymptotics we set out to derive 

P{x, t) = Pi°^ {ex - ewt, eh) + 0{e) . (2.28) 



3. Small-inertia expansion 

It is expedient to define the {d + l)-dimensional vector fields ip and W with compo- 
nent^ 

r — and Wf, := , for /i = 

I ■■= ^ and + I3u^{x, t) ~ , for 1 < /i < d . 

From (|2.15p and (|2.22qP it follows immediately that: 

£(°VM(a', V, t) = ~W^,{x, V, t)iboix, V, t) , (3.1a) 

V'^(a;,t>,0) = 5 f,op,j:{x)py{v) . (3.16) 

The problem of deriving systematic approximations for the auxiliary fields p'^^\ p*-^^ 
reduces to that of solving the spectral properties of L'^'^'^ in some suitable limit. The 
small-ine rtia limit recovering the over-damped Fokker-Pl anck from Kramers' dynamics 
(see, e.g., Gardiner 1985HRiskenlll989l : IVan Kampen 2007 ) is a natural candidate. To set 



the scene, let us denote by U the typical magnitude of the advecting velocity field {U 
may be defined e.g. as the space-time average of ||ti||(a;,t) over [0,T] x B). From [/, the 
unit cell linear size the Stokes time rs, the Brownian diffusivity k and the gravitational 
acceleration we can construct the hydrodynamical numbers 

TS.U U £U 

Stokes: St := ^— , Froudc: Fr := —= , Peclet: Pe := — . 

£ ' K 

f Notice that i/jq thus turns out to have dimensions (length"'^ x velocity"'*) different from 
those of the remaining components (length^~''xvelocity~''). One could easily fix this point by 
setting ipo ■= £p^'^\ but we prefer to avoid such an operation for the sake of notational simplicity. 
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The small-inertia limit corresponds to let St tend to zero for arbitrary but finite values 
of Fr and Pe. In such a limit the Stokes time sets the shortest time scale in the problem. 
The corresponding typical co-velocity variation is diffusion dominated. The root mean 
square along any Cartesian axis, 



StPe 



implies that co-velocities become fast variables in units of U for time scales of the order 
of Ts. We therefore attune the functional dependence on St to the description of slow 
degrees of freedom by working with rescaled co-velocities, 



St Pe 



(3.2) 



and by decoupling the time dependence of solutions of p. lap into a fast [t^) and a slow 
{t') component, so as to deal now with a new vector field ip*{x, v', t' , t^), which is trivially 
related to its original counterpart xp through the multiplication by the Jacobian. 
Upon inserting 

did d 
dt ~ StW ^ dt' 
into (|3.1ap . the variables (|3.2I) yield for the expansion 



jrio) ^ J2 St"/2-i£(0:«) 



n=0 



where (omitting from now on the prime over-script to simplify the notation); 

d 



£(0:0) .= _^ _ --—-V, 

dti Ti Ti dv ' 



(3.3a) 



£(0:: 




(3.3&) 



ot ox 



£(0:3) 



PeG 



d_ 
dv 



(3.3c) 



(3.3d) 



Here and in what follows, the digit after the colon in the superscripts represents the order 
of the expansion in St, keeping in mind that the latter is in half-integer powers and for 
some quantities might not start from the zeroth order; this notation was used in order 
to distinguish such an expansion from the multiscale one shown in the previous section. 
We defined in (|3.3ap 

2 ^ 

to identify the variance of the co- velocity steady-state Maxwell-Boltzmann distribution 
(see § lA.ll below) and the advective time scale, and in p.3(i[) 

iidii 



to exhibit the Froude number dependence. 
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Two remarks are here in order. First, the effect of the rescahng (|3.2p is to balance 
drift and diffusion terms in the co- velocity dynamics for vanishing St. The solution of 
the auxiliary problem is thus turned into a multiscale expansion in powers of St based 
on the well-known Ornstein-Uhlenbeck operator (13.3 ap (see § [B] for further details). 
The kernel of defines a Maxwell-Boltzmann co-velocity distribution to which any 

integrable initial density relaxes for time scales much larger than the Stokes time rg. 
Second, the functional dependence on the faster time is only of conceptual importance. 
It only describes the exponentially fast equilibration process of the marginal co-velocity 
distribution to the aforementioned Maxwell-Boltzmann steady state, and can therefore 
be neglected in the asymptotics we are interested in. The upshot is that we can restrict 
the attention to the t^-asymptotic form of the power-series expansion 



± 



n=0 



and take into account the initial condition (j3.1b|l only to the extent that it affects the 
time dependence upon the slower time variable. With this proviso, the expansion of the 
terminal velocity (|2.19p in powers of St yields 



(■_!) oo 

^ + VSt"/2«;(-) 



(3.4) 



where 



w 



i-.u) 




i-.n+l) 



i-.n) 



n> 



(3.5) 



so that the auxiliary equation (|2.22ap spawns a perturbative hierarchy the first terms 
whereof are (no sum on fi is implied here): 
0(St-^) : 



(3.6) 



0(St 



_£(O:l)^0O) 




^ (:-l) 



^0 



(:0) 



(3.7) 



0(St°) 



^(0:0)^(:2) ^_^^(0:2-„)^(:„) 



0) 



(3., 
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2 

£(0:0)^03) ^ _^£(O:3-n)^0») 



(3.9) 



We solved p.6p - (l3.9p by first deriving the coefficients of the expansion of i/jq = P^"^ 
(the phase-space density) to thereafter use them for the evaluation of the remaining 
components of the auxiliary vector field '0>i = p'-^\ The details of the calculation can 
be found in §|3 more precisely in ij lA.ll for the first step and in ij IA.2l for the second one. 

3.1. Expression of the eddy difjusivity 

Within third-order accuracy, the co- velocity dependence of both the phase-space density 
p^^^ and the auxiliary vector field p*^^^ is parity defined (see § In consequence, the 
Stokes-number expansion of the eddy diffusivity (|2.26p tensor reduces to: 



r 



dt 



-St 5/ ^ /da3 / Av\ W^z;^^(^3)-K/3u^^02) 

' JK JR<i \ V -re 



(3.10) 



Inserting the terms of the expansion of p^^-' leads to further simplifications (see § IA.3p . 

3.1.1. Evaluation of the 0(St*') term 

The tracer limit of the particle eddy diffusivity (first line of p.lOp ) reduces to 

= J^dxu,{x,t)e-"H=r^,t) , (3.11) 

with the auxiliary vector ^(i-'^' obeying an advection-diffusion equation with zero initial 
condition; 



d / X 5 . 

— +u[x,t) ■ kA^ 

at ox 



^'■'■■°\x,t) 



u{x, t) 



|(i^"Ha;,0) = 



An alternative rewrite of p. lip . 



r 



dt 



T 



dxe''§^ix,t)^^{x,t) 



(3.12a) 
(3.126) 

(3.13) 



shows that this quantity is positive definite and recovers the result of iBiferale et 
(|l995h . 

3.1.2. Evaluation of the 0{St^) term 

The leading-order correction in p. 101) turns out to be 



T 



dt 

r 



dx 



Ui^{x,t)^l^'-^\x,t) - '^^^^ ^ ^K f,{x,t)u^{x,t) 



(3.14) 
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where the auxiUary vector ^(i-^) satisfies another (more complicated) advection-diffusion 
equation with vanishing initial condition: 



<xj Ky ij-f 

9 ^(0:2) , d 



5 



Here, £("^2) 

was defined in (|3.3cp and 

1 -/3 , , [Yd 

M := \ —u{x,t) - aid — — 

T£ae V 2 V Pe oa; 



C(1:0) 

(3.15a) 
(3.156) 

(3.16) 



3.1.3. Perturbative expression of the eddy-diffusivity tensor 

Gleaning p.l3p and p.l4p , we get the final result for the Stokes- number expansion of 
the eddy-difFusivity tensor 1)2.26^ : 



Kfj_^ = k \ 5 



r 



^^tJC/ -y X 



r 



dt 



(1:2) 



+ 0(St2), (3.17) 



with ^(i-*^) and ^(i-^) given by p.l2ap and (|3.15a[) . respectively. An important physical 
implication of p.l7p is that, within 0(St)-accuracy, the eddy diffusivity depends only 
implicitly on gravity through p.lSap . 
A more compact rewriting is 

rT, 



K 



dt 

T 



dx 



-Su^£,l}^ + St(l - P)u^u, + 0{Se) 



with the auxiliary vector to be expanded in a power series in St, whose zeroth and 
first orders are ^(^-O) and ^f^-^'. 

3.2. Expansion at small inertia and finite terminal velocity 

The procedure exposed up to now assumes finite Fr, i.e. finite gravity eff ects, which in 
the limit of small St makes w infi nitesimal. An alternative point of view (jMaxevlll987 : 



Marchioli. Fantoni k, Soldatill2007l ) consists in assuming the terminal velocity as a finite 



parametre: we will investigate this in the present subsection, denoting with a tilde the 
quantities which undergo some change. 

A limit of physical interest is thus that of vanishing St for large values of Fr, so that 



rr 



u 



The limit implies a reshuffling of the hierarchy 
becomes same order as (j3.3bp . so that 



(3.18) 

since the operator ([33^ 



£(0:l) ^ ^(0:1) ,^ l^J^ [u{x^t) + S^UG] .^ + Mvf 

T£ \ 2 ov \ Pe ox 
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whilst 



^(0:3) ^ ^(0:3) ^ 



The redefinitions affect both the phase-space density p*^"-' and the auxihary field p*^^^ (see 
§ IA.4p . Plugging these results into the St-expansion of the eddy diffusivity, we obtain 



T 



dt 



Sf,u+ I y / da; 



5^(1:0) 



{x,t) 



(1:0) 



ix,t) 



+ 0(St), (3.19) 



which now exhibits an implicit dependence upon gravity already at O (St°) (for this 
reason here we do not write down the results for O (St^)) through the new following 
advection-diffusion equation, again supplemented with zero initial condition: 



- + [uix, t) + (1 - I3)SfUG] • ^ - '^A, 



iil-.Q) 



{x,t) 



u{x, t) 



(3.20a) 



(3.206) 



4. Specific examples: parallel flows and Kolmogorov flow 

A flow is defined parallel if it is always and everywhere in the same direction (e.g. along 
a;i), i.e. if the velocity field can be written as u^{x, t) = Sfj_iu{x, t). The incompressibility 
constraint of zero divergence then requires u to be independent from xi, i.e. 

Ufj,{x, t) = 5piiu{x2, ... ,Xd,t) . 

(Further simplifications can be obtained for specific instances of flow, e.g. if u^(x,t) = 
5^iu{xd,t), only dependent on one spatial direction along which we orient the Xd axis, 
as is the case for the well-known Kolmogorov flow.) 

Equations (|3.12ap and (|3.15ap . coupled with vanishing initial conditions for ^(i-") and 
make s}^''^\x,t) and s}i"'^\x,t) the only non-zero vectorial components, as they 
are the only ones forced. Moreover, this forcing is independent of xi, and thus the same 
independence also holds for S,^'*^^ {x2^ . . . , a^d, t) and ^J^'^' {x2, . . . , Xd, t). This means that 
the advective terms of the type u ■ d/dx always vanish, and that the aforementioned 
equations, which turn out to be heavily simplified, can easily be solved in Fourier space 
thanks to the absence of any convolution. We denote the Fourier transformed by a hat 
(for the sake of simplicity, the cell size is assumed as i in every spatial dimension), 

f{ng,n^) ^ r^r-^ l\t J^dxe-'^^-^+^'^fix^t) , 

with q := 2'Knq/£ and uj := 27rn(^/T, so that 

f{x,t)^ ^ ^ e'('^-+"*)/(n„nJ. 

Therefore, we have: 

[at ^^y^^ - ^rf"=^^M ("q,"-)- gd i^ + ^g2 ' (4.1) 
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dt 



t(l:2) ^ 



d 



d 



K-'TfAi - i^n-g^A^ -Te{l- (3)g ■ — 



dx 



.(1:0) _ n{l-(3)K 



(1 + P)iTeK'^q^ + iTiKUjq^) - 2iTf (1 -p)g-q 



2{iuj + Kq^Y 

Upon antitransforming to the physical space and plugging into p.l3p and (|3.14p . we get: 



Jo I JK 



X 



p ■ q 



,J[{p+q)-x + (i-l+ui)t] 



(4.2) 



Jo I Jb l „ „ „ 



(1 + /3)(T^«^g^ + ir^At^g') - 2iT^(l - /3)g ■ g 



+q)-x+(J2+l^)t] 



2(ia; + Kq"^ 



T 



dt 



T 



da; T£(l — P)u^{x, t) 

, ,2 (1 + ;3)(r^K294 + innLoq^) - 2ir<>(l - fi)g ■ q 



2{\UJ + Kg2)2 



(4.3) 



For generic parallel flows, gravity thus does appear in the 0(St) correction to the eddy 
difFusivity, (|4.3p . and can in principle bring about either an enhancement or a depletion of 
the tensorial component along the flow directio n. The tracer hmit of the eddy difFusivity, 
(|4.2p . obviously coinc ides with the results of Biferale et al. ( 19951 ) . and always shows 
a positive correction ( Pavliotis fc Stuart 20051 ) with respect to the Brownian isotropic 
counterpart in the absence of any flow. 

4.0.1. Steady Kolmogorov flow 

These result s can be s ubstantiated even more deeply by focusing on the steady Kol- 
mogorov flow (|Obukhovl[l98a) 

Uf,{x, t) = Sf,iU cos{qeXd) , (4.4) 
with qi := 2'k/£, such that u{nq,n^) = C/(5„„^o<5„,^ ,o ■ ■■5n,^_^,o[5n,^,i+5n,^-i)/'i- Thus, 



C/2 

K^u = + — ^(5^i(5,yi + St 
^q^ K 



In other words. 



i^ll := Kii = K 



Pe" 



t/£(3 - /3) 



3-/3 



St^iSui + 0(St2) 



1 + TT^ + — r^PeSt + 0(St' 



(4.5) 



(4.6) 



i.e. the component of the eddy diffusivity parallel to the flow is enhanced by the presence 
of the latter, and gravity can play a role only in subleading corrections at higher orders 
in St, independently of the cell orientation. This enhancement is maximum at larger and 
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larger Pe already at the tracer level, and for very heavy particles taking into account 
the role of inertia at its leading order as well. The further enhancement due to inertia 
can easily be explained by reminding that, according to the Taylor formula, the eddy 
diffusivity is loosely proportional to the particle Lagrangian autocorrelation time, which 
is definitely longer for inertial particles than for tracers because of ballistic "coherence" 
or "memory" (notice that this mechanism provides an explanation for the present steady 
flow, but in general need not hold in unsteady cases). On the other hand, no enhancement 
is present for a vanishing Peclet number, nor for bubbles. 
The perpendicular component is not modified within accuracy, 

^^j/"'" ^.[l + 0(St^)] , (4.7) 

and it is worth noticing that the Xd axis does not introduce any preferential direction in 
the {d — l)-dimensional subspace orthogonal to the flow, in which the eddy diffusivity is 
isotropic. 

4.0.2. Random-in-time Kolmogorov flow 

As a further example, let us consider a random-in-time Kolmogorov fiow, 

u^{x,t) = Sf,iU cos{qeXd)F{t) , (4.8) 

such that the random te mporal part is statistically stationary with correlation function 
( Monin fc Yaglon]|ll975l ) 

C(iA) {Fit + tA)Fit)) = e-"l*^l cos(7tA) • (4.9) 

This expression requires the introduction of two additional adimensional numbers, 

Kubo: Ku :— , Strouhal: Sr ^^-^ , 

are ztt 

which give a measure of the vortex life time and of the recirculation degree, respectively, 
in advective-time units (notice that, in the spirit of a comparison with the time-periodic 
case, 7 can roughly be identified with 27r/T). The formalism we applied up to now still 
holds, with the proviso of dealing with a transform rather than a series in the Fourier 
frequency domain, and of replacing any time average over a period with a statistical 
average (•). We have indeed 

— ^ M(np,no)u(nq,n<^)e''"+'^'*/('^:^^) / da;|M*(nq,w)|V('^,-w) 



with |uH.(nq,w)p = U'^F{uj)Sng^fi ■ ■ ■ 5n„^_^,(){Sng^,i +^ri,^ -i)/4 and 

1 1 



Pico). " 



a2 + (cj + 7)2 + {oj - 7)2 

From (|4.2p - (l4.3p . the upshot is a modification only of the parallel eddy diffusivity: 



" \\ 87r2 Pe 



1-13 l + (3Pj 



K 



PeSt-hO(St2) } . (4.10) 



Here, 

p _ 4^2 a + qj^ _ Pe-^+Ku'V4^ 

' n{a + q}nf+-i^ (Pe"^ + Ku" V4^2)2 + SrV47r2 ^' ' 



16 



M. Martins Afonso, A. Mazzino and P. Muratore-Ginanneschi 



1 




0.8 




\ 
\ 
\ 


0.6 












0.4 






0.2 








Ku = 


l = Sr '\ \ 




Ku - 


100, Sr = 0.01 , \ 




Ku = 


0.1, Sr - 10 " • 








0.01 0.1 1 10 100 1000 10000 

Pe 



Figure 1. Ratio Px/Pe as a function of Pe for different values of Ku and Sr. Notice that the 
ratio is closer and closer to 1 for large Ku and small Sr, in which case 1)4. 10|) tends to (|4.6p . 



can be identified as a "modified" Peclet number, which gets back its "original" value 
Pe in the steady case a — >■ 7, i.e. for infinite Ku and vanishing Sr (in which 
case (|4.10|) coincides with (|4.6|) ). In general, a non- vanishing decaying rate a (finite 
Ku) helps increasing the role of the molecular diffusivity in the expression of the eddy 
diffusivity, while a finite recirculation degree 7 (non- vanishing Sr) contributes to reducing 
the parametre Pk- 

The outcome is that, at the tracer level (St — )■ in (|4.10p ). the correction with respect to 
the Brownian value is always positive, but its magnitude is smaller with respect to the 
steady-Kolmogorov correction (|4.6p because of a multiplicative factor Pk/Pg (< 1 and 
investigated in figures [THU) . 

For what concerns the role of inertia at its leading order, again the magnitude of the 
correction is smaller with respect to the steady-Kolmogorov case (note that (3 — /3)/4 = 
(1 — l3)/2 + {1 + 13) /4:, for a comparison with ()4.6p ). but its sign turns out to show a 
diffusivity depletion for particles light enough, as the coefficient in square brackets in 
(|4.10p is negative for (3 > (3* . Here, the critical value 



2 - P^/Pe 



e [1.5,3] 



is a growing function of Pk, and thus gets back its "steady" value 3 for infinite Ku and 
vanishing Sr. 



An interesting quant ity to investigate is the interference (jMazzino fc Vergassolalll997 : 
Castiglione et al. 19981) between the different physical mechanisms into play, which pos- 
sess different intrinsic time scales (advective, diffusive, inertial, 1/a, I/7) and can thus 
interfere constructively — to give a resonance — or destructively. We then compute the 
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0.001 0.01 01 1 10 0.1 1 10 100 1000 

Kii Si 



Figure 2. Left: ratio Pk /Pe as a function of Ku for different values of Pe and Sr. Right: ratio 
Pk /Pe as a function of Sr for different values of Pe and Ku. 



St < 3.05 10""* St = 10"^-^ St = 10"^ St > 1.22 10"^ 



$13 13 > 2.86 /? > 0.22 V/3 

(no type) (only bubbles) (also tracers) (every particle) 

Table 1. Focusing on the case Pe = Ku = Sr = 1, i.e. a situation where the interference is 
destructive in the tracer approximation St — >■ 0, the table reports (for different, small values of 
St) the ranges of /3 which turn the interference constructive upon taking inertia into account at 
its leading order. 



quantity 



K.:= Kn -K- lim K« = U£ 



^(PK-Pl^) + ^^St + 0(St^ 



(4.12) 



where 



P* lim Pk = 5— ^ . (4.13) 

A positive or negative sign of K. then indicates a constructive or destructive interference 
of the Brownian diffusivity vifith the other mechanisms. 

At the tracer level, this corresponds to investigate the sign of (Pr- — -P^), which is done 
in figure |3] and shows the possible presence of both alternatives. Analytically, one finds 
a critical value of the Peclet number, 



Pe* 



Ku^i 



Sr^ - Ku-V47r2 ' 



such that the interference is constructive for Pe > Pe* and destructive otherwise (the 
interference is always destructive if Pe* turns out to be negative, e.g. for vanishing Sr). 
Then, one can also take inertia into account at its leading order, which always carries a 
positive addend in /C, whose magnitude increases linearly with both /3 and St. Starting 
from situations where /C < in the tracer approximation St —> 0, it is thus interesting 
to investigate for what ranges of /3 and St the interference becomes positive: an example 
is sketched in table [TJ 
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11, 



0.5 



-0.5 



0.001 
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Ku = 10, Sr = 1 
Ku= 1, Sr = 10 




1000 



Figure 3. Sign of the interference (|4.12|l as a function of Pe for different values of Ku and Sr. 
Note that the two dashed lines correspond to the same value of the product KuSr but give 
different results, which means that the ratio y/a from (|4.9p alone in this case is not sufficient 
to determine the nature of the interference. 



4.1. Expansion at small inertia and finite terminal velocity 

Let us now investigate the case of finite terminal velocity p.lSp in the framework of 
parallel flows, i.e. let us consider p.20ap instead of (|3.12a|) . We then obtain a modified 
version of (14.11): 



(1:0) 



i'i ILJ + 1(1 - P)SfUG -q + nq^ 
From p.l9p . we deduce that (|4.2p is to be modified into: 



/o TJm 



{p ■ q)e'[(P+'?)-^+(f*+")*l 



[in + i(l - P)SfUG ■ p + Kp2][icj + i(l - i3)SfUG ■ q + nq^] 

fLj + n - 



, c c V- \u{nq,n^)\^q^ 

+ d.id.i 2^ 2^ ^ _ ^^^^^^ ^ ^ 



4.1.1. Steady Kolmogorov flow 

Let us now impose once again the steady Kolmogorov flow (|4.4p and, instead of (|4.5p . 
we get 

[/2 



2f/2(l - /3)2S'2 cos2 e + 2qjK^ 



0(St) 
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where 9 is the angle between gravity and the Xd direction. In other words, (|4.6|) - (|4.7p 
now become 



87r2 



2Pe2(l-/3)25|,cos2 6' 



0(St) 



(4.14) 



K_i_^ K_i_ = k[1 + 0(St)] 



(4.15) 



Comparing (I4.14p with the corresponding result for vanishing terminal velocity, (j4.6|) . we 
deduce that the effect of gravity at 0(St*') is to reduce the enhancement of the parallel 
eddy diffusivity caused by the presence of the flow with respect to the Brownian value. 
This reduction of the eddy-diffusivity enhancement is most pronounced when the Xd 
direction (on which the Kolmogorov flow depends sinusoidally) is vertically aligned — 
independently whether upwards or downwards — , which also means that no enhancement 
reduction is present if the flow is vertical {xi \\ G ^ Xd -L G =^ cos 6 = 0). For what 
concerns the eddy diffusivity perpendicular to the flow, a comparison between (j4.7|) and 
(|4.15p shows no change at 0(St°), but a finite terminal velocity can in general move the 
effect of gravity from the original O(St^) down to O(St^). 

Alternatively, one may study the eddy diffusivity in the vertical and any horizontal 
direction; their expressions coincide with ()4.14|) but with the fraction in square brackets 
multiplied by cos^ 8 and {d— 1)~^ sin^ 8 respectively, where 8 is the angle between the 
vertical and the flow directions. 



4.1.2. Random-in-time Kolmogorov flow 

If, on the contrary, the random Kolmogorov flow ()4.8|) - (|4.9[) is imposed, the eddy 
diffusivity is again modified only in its parallel component, according not to (j4.10p but 
rather to 

i^ii i-> i^ii = K 

Here, (|4.1ip becomes 



a + qjn 
2 ' 

with 

r± := 7 ± Uqi{l - /3)Sf cos 9 . 

This latter expression shows that sedimentation contributes to change the flow character- 
istic frequency 7 as viewed (or "sampled") by the particles, but in a manner symmetric 
for up/down reflections (0 n- tt — 0) because of the symmetrization present in Pk (clearly 
no modification occurs for neutral particles /3 = 1, nor for the orthogonal case 9 ~ 7r/2). 
In the steady case a = = 7, expression (|4.16p obviously coincides with (|4.14|) . 
Another interesting comparison for (|4.16|) is with the corresponding result (|4.10[) for van- 
ishing terminal velocity (limit Sf — ^ 0); then the ratio Pk/Pk has to be investigated, 
which is done in figures HHSl the upshot is that, in the framework of the random-in-time 
Kolmogorov flow, gravity can either increase or decrease the enhancement of the parallel 
eddy diffusivity, according to the values of the other parametres. 

Let us now analyse the problem of interference again. From the usual point of view. 



(4.16) 



Pk ^ Pk 



47r 



2 



(a -f (7«k)2 



(a + qfny 
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Ku= Sr = 5j.^ = 1, y3 = = 9 
- Ku = 10, Sr = 1 = /3 = = 9 
Ku = 1 = Sp, Sr = 0.1, = O = e 
Ku = 1 = Sr, Sf = 10, ^ = = e 
Ku = Sr = Sf = 1, 3 = 3, 8 = 
Ku = Sr = Sf = 1. 3 = 0. e = tt/I 




- Pe = Sr = Sf = 1, ^ = = « 

- Pe= 10, Sr = 1 = Sf, /3 = = 8 
Pe = 1 = Sf, Sr = 10, ^ = = e 

- Pe = 1 = Sr, Sf = 10, 3 = = 8 
Pe = Sr = Sf = 1, ,3 = 3, 8 = 
Pe = Sr = Sf = 1, ,3 = 0, 8 = 7r/4 



Figure 4. Left: ratio Pk/Pk as a function of Pe for different values of Ku, Sr, Sf, P and 9; the 
ratio is closer and closer to 1 for small Sf, (1 — P) or {n/2 — 9), in which case (|4.16p tends to 
1)4. 10|) . Right: ratio Pk/Pk as a function of Ku for different values of Pe, Sr, Sf, P and 9. 




/ Pe = Ku= Sf = 1. 3 = = 6 
/Pe = 10. Ku = 1 = Sf, y3 = = 
•' Pe = 10 = Ku, Sf = 1, ^ = = 
Pe = 1 = Ku, Sf = 10, ^ = = 
Pe = Ku = Sf = 1. 3 = 3. 8 = - 
Pe = Ku = .Sf = 1. 3 = 0. 8 = jr/4 



100 









Pe = Ku = Sr = 1. 8 ={] = e 

Pe = 10, Ku = 1 = ar. 3 = = 

Pe = 1 = Sr. Ku = 0.1. /3 = = 

Pe = 1 = Ku. Sr = 10. /5 = = 

Pe = Ku = Sr = 1. ^ = 3. 6 = 
Pe = Ku = Sr = 1. 8 = 0. 9 = 


» \ 

\ 





Figure 5. Left: ratio Pk/Pk as a function of Sr for different values of Pe, Ku, Sf, P and 9. 
Right: ratio Pk/Pk as a function of Sf for different values of Pe, Ku, Sr, P and 9. 



Pe = Ku = Sr = Sf = 1. * = 

Pe = 10, Ku = Sr = Sf = 1, 8 = 

Pe = 10 = Ku, Sr = 1 = Sf, 8 = 

Pe = Ku = Sf = 1, Sr = 10, 8 = 

- Pe = Ku = Sr = 1. Sf = 10, 8 = 
,- Pe = Ku = Sr = Sf = 1, 8 = ir/4 . - 



- Pe = Ku = Sr = Sf = 1, ^ = 

- Pe = 10, Ku = Sr = Sf = 1, (3 = 
Pe = 10 = Ku, Sr = 1 = Sf, ;3 = 
Pe = 10 = Sr, Ku = 1 = Sf, /3 = - 
Pe = Ku = Sr = 1, Sf = 10. /3 = 
Pe = Ku = Sr = Sf = 1, 3 = 3 



Figure 6. Left: ratio Pk/Pk as a function of P for different values of Pe, Ku, Sr, Sf and 
Right: ratio Pk/Pk as a function of 9 for different values of Pe, Ku, Sr, Sf and p. 



this implies investigating not ()4.12|) but rather 



IC ^ K. := Kn - K ~ lim Kn = U£ 

" Pe-J-oo " 



87r2 



(Pk - P*k) + 0(St) 



(4.17) 
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Ku = Sr = 5"^ = 1, /3 = 0, 6 = 7r/4 
Ku = Sr = 5f = 1, /3 = 3, 6> = 
Ku = 1 = Sr, 5j. = 10, I3 = Q = B 
Ku = 1 = Sf, Sr = 0.1, /3 = = e 
Ku = 0.1, Sr = 1 = 5^., /? = = 9,i 




1000 



Figure 7. Sign of the interference (|4.17|) as a function of Pe for different values of Ku, Sr, Sf, 

13 and 9. 



with (I4.13P now becoming 
Pk^Pk-^ lim Pk 



Po— )-CX) 



Ku 



-1 



2 1 Ku" V47r2 + [Sr + (1 - P)Sf cos 6l]2 Ku" V47r2 + [Sr - (1 - P)Sf cos 6] 

A positive or negative sign of K, then indicates a constructive or destructive interference 
of the Brownian diffusivity with the other mechanisms: this corresponds to investigating 
the sign of [Pk — Pk) as a function of Pe (oc which is done in figure [7] and shows 

the possible presence of both alternatives. 

However, one could also wonder how gravity interferes with the other mechanisms. There- 
fore, one can study the sign of the quantity 



^l|l.=o 



1=0 



lim 

Pc— >o< 



(^11 -^ll) 



Ul 



1 

8^ 



.{Pk 



P, 



K) 



0(St) 



(4.18) 

i.e. of (P^ — P^) as a function of Sp oc ||g|| (where the Brownian diffusivity must be 
neglected for sake of consistency). This is done in figure H] and again shows the possibility 
of having either destructive or constructive interference, with the relevant annotation that 
the latter can now be very pronounced near well-defined resonance peaks at critical values 



S*F 



|1 -/3|cos( 



(except for small values of Ku, which make the denominators in P^ large anyway) . 
Notice that the situat ion analysed in t he present subsection i s a gen eralization of the 



iNotice tnat trie situat ion analysed m t ne present subsection i s a gen eralization ol tne 
one already attacked in 'Mazzinol ( 19971 ): iMazzino fc Vergassolal (1997), with the role of 
the mean, large-scale drift played here by sedimentation. 
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Figure 8. Sign of the interference (|4.18|l as a function of Sf for different values of Ku, Sr, /3 

and 6. 

5. Conclusions and perspectives 

We have investigated the asymptotic (in space and time) behaviour of the concentration 
field of inertial particles of arbitrary mass density carried by an incompressible flow field 
(static or time dependent) under the action of gravity. 

Exploiting the multiple-scale perturbative expansion in the scale-separation parametre, 
we found a diffusive equation for the asymptotic dynamics of the particle concentration 
field, involving a positive definite eddy-diffusivity tensor. The method reduces the com- 
putation of this tensor to the solution of an auxiliary differential problem valid for any 
Stokes, Froude and Peclet numbers. The general results have been greatly simplified 
focusing on two relevant limits of vanishing inertia: one with vanishing sedimenting ve- 
locity, the other maintaining a finite value of the latter quantity. 

As also happening in the simpler problem of inertialess-particle transport, the deter- 
mination of the eddy-diffusivity tensor is in general possible only exploiting numerical 
strategies. Here the problem is particularly evident owing to the larger (compared to 
the inertialess case) phase space. Nevertheless, we were able to select a relevant class of 
carrier fiows (parallel flows, both static and time-dependent) where, after some manip- 
ulations, one arrives to exact, explicit expressions for the eddy diffusivities. The latter 
expressions allow us to extricate a highly nontrivial role played by gravity and inertia on 
the diffusive behaviour, where regimes of both enhanced and reduced diffusion can be 
observed depending on the detailed structure of the flow. 

Here is a summary of our findings for the parallel Kolmogorov flows, highlighting the 
nontrivial role of inertia and gravity on the diffusion process in the direction parallel to 
the flow. 

• An enhancement of the diffusivity is always found, and the Peclet number always 
plays a positive role in this respect (more precisely, analysing the ratio K\\ /k). However, if 
one compares the situations in which other mechanisms — such as flow time dependence 
or sedimentation or both — are present, with the corresponding situations in which the 




Ku = 1 = Sr, l3 = Q = e 
Ku = 1 = Sr, /? = 0, e = 7r/4 
Ku = 1 = Sr, /? = 3, e = 
Ku = 1, Sr = 10, /3 = = 61 
Ku = 0.1, Sr = l, 13 = = 6 
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latter are absent, one can find an enhancement reduction which is more or less pronounced 
at different Pe. By introducing the concept of interference, one concludes that Brownian 
diffusivity interferes constructively /destructively with the other mechanisms for Peclet 
numbers larger/smaller than a critical value, respectively, but the latter can in principle 
be made arbitrarily small or large. These conclusions hold in general, irrespectively of 
considering a tracer or an inertial particle at small inertia, endowed with vanishing or 
finite sedimentation. 

• The role of inertia (Stokes number) should be taken into account together with the 
coefHcient /3 based on the density ratio. The diffusivity enhancement indeed shows a 
leading correction proportional to St, which is always positive (and decreasing with /?) 
for a steady flow, but which can become negative for a random time-dependent flow if the 
particles are light enough. This enhancement reduction, while passing from the former 
kind of flow to the latter, is a decreasing function of the Strouhal number — i.e. vortex 
frequency or recirculation degree — and a growing function (loosely speaking, except for 
some peculiar cases or ranges) of the Kubo number — i.e. the life time of the structures 
before decay. (Also notice that Ku always appears in the form of a harmonic sum with 
Pe, as resistances in parallel.) St and /3 also contribute linearly and positively to the 
above-mentioned interference. 

• Gravity (Froude number) plays a significant role only through sedimentation, i.e. 
when the particle terminal velocity {Sp, modulo a prefactor) remains finite even for 
vanishing inertia. In that case, the key quantity is the combination (1 — I3)SfCos9, 9 
being the angle between the vertical and the Xd direction upon which the flow depends 
sinusoidally (and NOT the xi direction along which the flow points: this distinction is 
crucial in the three-dimensional case) . This key parametre brings about a monotonically 
growing reduction of the diffusivity enhancement in the steady case. In the random time- 
dependent case, this enhancement reduction: coincides with the vanishing-sedimentation 
value for a vanishing deviation of /3 from 1, or of S'i? from 0, or of from n/2; is maximum 
for some intermediate deviation, corresponding to a resonance with the flow frequency Sr; 
and vanishes for very large deviations, i.e. infinite Sp- Such a behaviour is summarised 
by the presence of strong resonance peaks denoting constructive interference of gravity 
with the other mechanisms; this interference becomes destructive for Sp large enough. 

We thank Antonio Celani for useful discussions and suggestions. MMA acknowledges 
financial support from "Institut de Radioprotection et de Siirete Nucleaire" , and PMG 
from "Center of Excellence AnDy of the Finnish Academy" . 



Appendix A. Details of the resolution of (13. 6p (13.91) for 

A.l. Phase-space density -00 = P^^^ 
The Stokes-number ex pansion of the phase-space density was studied in details in 
Martins Afo"^ (|2008l ). the main results whereof we hereby shortly revisit by couching 



them into our notation. The results the Ornstein-Uhlenbeck process recalled in §|B]imply 
that any integrable initial co-velocity distribution evolving in t^ according to (I3.3ap 
converges exponentially to the Maxwell-Boltzmann distribution 

pi^-^\x,v,t) = £^1^^(0:0) ^^^^^ ^ /-(«)e(°^°)(a.,0 . (Al) 

The function ^(o^") is a constant with respect to the fast variables {v,t^). We leverage it to 
impose the solvability conditions of the perturbative hierarchy starting at ()A ip . Higher- 
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order corrections are most conveniently evaluated observing that the grading (|B 5|) and 
recursion ()B 6|) relations, enjoyed by the Hermite polynomials i/"'°'*'s and their L^(K'^)- 
orthogonal duals x"''^*'s, imply that for any f{x) 

k=l 



oik 



an 



C'-'-'hZ.to.Mfix) = -J^ (L^l^;,«+i;^-^ (^)/(^) . (A3) 

We used for the x's the conventions fixed in §|B1 and we defined Ad through (|3.16p . 
The leading-order correction does not impose any solvability condition. Since by (|B 7p 

acts on ^"■'•'^^ as a multiplication by Ti/n (asymptotically in t^), we find 

p^°--^Hx,v,t) = r.xif (t;)X„,C'"^")(^,0 ■ (A 4) 

The function ^(0'°) in (jA II) and (|A4[) is specified by the solvability condition brought 
about by p.Sp . the solution whereof is 



{x, V, t) = x"^^' {v)£.^'--^^ {x, t) + ^xia (^)X.,A^„,^("^") (a., i) (A 5) 

if and only if ^(''■'') satisfies the advection-diffusion equation 

AAe'"^")(a;,t) =0, (A 6a) 

^("^")(a;,0) = p(a;), (A 66) 

with 

. r , . d . 5 , , 9 2rf (7? , 

"^■■-dl^^^^^'^-d-x-^^-^dl^^^^^'^-d^-^^- 
The solvability condition (jA 6aP is the Fokker-Planck limit for incompressible advection 
of the Kramers dynamic s for vanishing Stokes. As usual in the non-periodic case 
( Pavliotis &: Stuart 20051 ). the solvability condition arises from the requirement that the 



non-homogeneous term in (|3.8I) be orthogonal to the kernel of the zero-frequency compo- 
nent of £(0^°)^ i.e. the generator of the Ornstein-Uhlenbeck process ()B with respect to 
the scalar product on L^(]R'') (equation (jB SP ). The fulfillment of this condition prevents 
the occurrence of marginal (polynomial in t^) instabilities in the perturbative expansion. 
By (jA 6al) , any integrable initial condition on the periodicity cell B converges to a uni- 
form distribution. Since we are ultimately interested in degrees of freedom slower than 
the t-dynamics, it is not restrictive to choose as initial condition the invariant measure 

e("^")(a;,t)=p.(a;) = l. (A 7) 

Finally, the presence in (|A b\ of the homogeneous term x^'*^'^*^"'^' takes into account the 
solvability conditions intervening in the equation for p^*^-**) (similar to p.9p with each 
superscript to tp augmented by one). Namely, parity considerations on the co- velocity 
distribution evince that p.9p in its original form does not require any solvability condition 



f throughout the article, unless explicitly stated, the Einstein convention of implicit summa- 
tion holds over repeated indices Oi = 1, . . . , d (i = 1, . . . , n, for any n). 
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and yields 



P 



(0:3) 



^.'(1-/3) 



having used (|A 7|) and 



'2 3a/ 



1-/3 /Pe 



TiCFlt 



Ua{x, t) 



ait"- V 2 
d 



(0:2)Mqi 



The complete expression of p'*') within third-order accuracy in St is then attained by 
applying (|A 2p to (lA 8t, which shows that ^(O'^) satisfies the following advection-diffusion 
equation in B fjMartins Afonso 20Q8t ): 



a 1 ^^02 



(A9a) 



e(°^')(a;,0) =0 . (A9&) 
Owing to periodicity and incompressibility, the forcing sustaining (jA 9al) is zero average, 



dx. 



-{x,t)^^{x,t) = - I dxua2{x,t) ^ {x,t) 



= 0, 



(A 10) 



SO ensuring probability conservation within the consistent accuracy: 

/ dxC^°--^Hx.t) = / ^ /da; / dvp'^'>--^Hx,v,t) = . 
Jm Jo ' Jm Jr'' 

A. 2. Auxiliary vector field ip^ = p^jl\ /i = 1, . . . , d 

The Stokes- number expansion of the auxiliary vector field p'^' differs from that of the 
phase-space density p^^^ only by the presence of extra non-homogeneous terms. These 
latter simplify when the coefficients p.Sp of the Stokes expansion of the terminal velocity 
are evaluated using the results of the previous subsection for the phase-space density. In 
particular we have 

i^(-i) = = , 

since the co- velocity dependence of the coefficients of the phase-space density expansion 
is defined under parity 

p^^--''\x,-v,t) = (-l)V"^")(a;,i;,t) 
and by (|B4)) . By ([^ we also have 



(All) 



w 



T 



(-) = / ^ / dW d. ^ ^il-^^^iv)u.^ix.i) ^ ^^(^^ (.) 



dt 



r 



Id 



= 



The first non-vanishing contribution is then 

rT 



,(:2) 



f ^ fdx u{x, t)^^"--^^ {x, t) + li)g , 



using (|B 4p again and the vanishing of integrals (over the period) of derivatives of periodic 
functions. 
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A. 2.1. Leading orders of the Stokes expansion 

Proceeding as for the phase-space density, we look for a solution of (|3.6I) in the form 

pl'--°\x,v,t) ^ X"'''''{v)^l'--"\x,t) . (A 12) 

In full analogy to § lA.l) the vector field ^(i-*^) is determined by the first solvability 
condition brought about by the perturbative hierarchy. The analogy with the phase- 
space density extends inasmuch that remains bounded when acting on the 
non-homogeneous terms sustaining the first correction term in Stokes (j3.7p . whence 



p';!--'Hx,v,t) = nx];^'{v)M^,^l'--"Hx,t) 



(A 13) 



V Pe ■ 

Thus, is specified by the solvability condition associated to the second-order equa- 

tion (|3.8p . This latter differs from (jA 6a| by the occurrence of a non- homogeneous term 
on the right-hand side, and gives p.l2ap complemented by the initial condition (j3.12&p . 
An important consequence of the very definition of terminal velocity (j2.19p preserved by 
its Stokes-number expansion p.4p is the conservation law 

[dx^^^-°\x,t)= [dx [ dvp''^--^^\x,v,t) =0 . (A14) 

Jb Jo I Jm Jr'' 



(A 15) 



Inserting (|3.12a|) into (|3.8p yields 

pi'--'\x,v,t) = x'-''''{v)^i'--'Hx,t) 



-y (v) 



Like (jA 5|) . (|A ISp comprises a term x^''^'^ ^(i-^) solution of the associated homogeneous 
equation. The introduction of such term is required to impose the solvability condition 
in the equation for p^^-^^ (in which, with respect to p.9p . beyond the aforementioned 
augmentation of the superscripts, an additional term — iw^'^'^T/'g °' appears within square 
brackets) . Again, parity considerations (|A lip as in § lA.ll entail that equation p.9p does 
not engender any solvability condition. Tedious but straightforward calculations using 
((X2)) . (|X3l) and (HH) yield 



Fr^ 



2 l;<Tf/'(0:2)A^ t(l:0) , 3;fTf 



(1:0) _ 



r|(l-/?)_ 



T£/3(l - /3) 



Id 



dx,, 



The expression within third-order accuracy of p'^^-' is complete observing that ^(i-^) sat- 
isfies the solvability condition (I3.15al) with initial condition (|3.156p . and the conservation 
law 

/ dxi^^--^Hx,t) ^ [ % Idx ( dvp'^^-^Xx.v^t) = . 
Jb Jo ' Jb Jh^' 
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A. 3. Expression of the eddy difjusivity 

A. 3.1. Evaluation of the 0(St°) term 

Using (lA 121) . (lA 131) and expressing the co- velocity dependence in the integrands as 
the L^(IR'')-scalar product between the Hermite polynomials iJ"''^* and their duals 
we have 









dv 


Jo ' Jk J 



T 



dt 



= -Sj ^ Idx I dvHl^'^'iv) 



(1:0) 



-sT^ [dx [ dvH°-''^^iv)(3x°--'''{v)u,e'-°^ 

Jo I JM JRt' 



The orthogonahty relations (jB 4[) . together with 



r 



y / da; 



(1:0) 



yield p. lip . The positive definiteness of Kii^ can be proven through p.l2ap and 



(:0) 



-{x,t) = 



r 



dt 



r 



da; 



A. 3. 2. Evaluation of the O(St^) term 
The leading-order correction 

r'^dt 

lo T 



dx / dv 



can be evaluated as the sum of two distinct integrals. 

The first involves the projection of p'l}'^^ over Hj^''^". Using (|B 4p to integrate out the 
co-velocity dependence, we are left with 



f %!dx[ dvaJ^H];-^{v)p^^-^) 
Jo I Jr Jr-^ V i^e 



•^dt 

r 

"S/da; 



T 



Tjaj_ _o_ 
Pe dxa 

2 



Pc 2£'' 



T 



dt 



-Sj^^ldx 



(A 16) 



where the latter simplifications are due to the conservation laws (jA 10| and (jA 14^ , com- 
bined with the vanishing of integrals (over the period) of derivatives of periodic functions. 
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The second term involves (jA Isp and is more straightforward to evaluate: 

Jo I JB JR-i Jo ' Jb ^ ^ 

(A 17) 

Expressions (jA 16[) and (jA 17|) combine into (I3.14p . which can be gleaned with p. 131) to 
give the final result (j3.17l) . 

A. 4. Expansion at small inertia and finite terminal velocity 
The leading-order correction in Stokes to the phase-space density becomes 

p(0:i) ^p(0:i)(^^„^i) = nx];:'{v)M^,i^"--"\x,t) , (A18) 

where now ^(0'°) satisfies 

Xi^°'-°Hx,t) = 0, (A 19) 

with 

d d d d Ue 

The new solvability condition (jA 19p remains however compatible with the steady-state 
solution (jA 7p . which we therefore assume. A further consequence of IjA 18j) is that the 
first non- vanishing contribution to the terminal velocity occurs at 0(St°): 

U;(:0) ^ ^{-.0) = (1 _ P)SfUG . 

The leading-order correction to p^^^ on its turn becomes 

^ p(^--^)[x,V,t) = nx];:^{v)M^A('--'\x,t) - I^^^X'-'^^ , 

where now ^(i-"' obeys p.20ap supplemented with (j3.205l) . Correspondingly, throughout 
this section (13. 16^ becomes 

(1 ~ B) fPe I 2 d 

M^M:=^- —\hr [uix,t) + SFUG]- J—aeTT ■ 

Tiui V 2 V Pe aa; 



Appendix B. Ornstein Uhlenbeck process 

We briefly recall some properties of the well-known Ornstein-Uhlenbeck process (see, 
e.g.. lVan Kampenll2007l and references therein) used in the main text. The generator of 
the Ornstein-Uhlenbeck process on with parametres r, cr > 0, defined as 

„2 



Ql. 



-A, 



y 

T dy T 

is diagonal on the basis of the multidimensional Hermite polynomials iJ^j' 



(Bl) 



^ T,a ai . ..an 



_ TTn;a 
. a 1 ... a ,5 



with the 1 < tti < d (i = 1, . . . , n) indices ranging over the y-vector components. Hermite 
(multi-)polynomials are most conveniently defined by the Rodrigucs formula 



Kr...„(y) = (-ir 



n 



d 
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where 



(27r(T2)rf/2 

(y2 = ||y|P) is the unique element of the kernel (steady-state measure) of the adjoint 
gr,a operator on L2(M''): 



yT,a ^ ^ h — Ay = X yT.a—„ 

The last equality in (|B 2[) implies that Qr.a is diagonal, 

n 



(B2) 



VT,aXai. 



'^ai...an ' 



on the L2(R'') complete set 



We also use the convention 



(B3) 



X2r.a„=^^"r..a„=0, ifn<0. 

The Rodrigues formula implies the L2(M'')-orthogonality relations between Hermite poly- 
nomials and the basis elements (jB 31) . Namely, a multiple integration by parts readily 
shows that the integral is non- vanishing only for ni — 112 — n, and if there exists at least 
one permutation V of the n-tuple . . . , such that ai — Vi{Pi, . . . , Pn)' 



d 



^ ( fi ^1 



— ^mna X/ n '^ai,-Pi(/3i,...,/3„ 



(B4) 



(the sum ranging over all permutations of the (/3i . . . /3n)). Further consequences are the 
grading relation 



^Aqi...q„ Aqi 



and the recursion relation 



(B5) 



(B6) 



which follows by applying (jB 4[) and (|B 5|) to 



Finally we observe that 



TTn-\-l-(T (li'] 
"(y) ~ 1" o J^a[...a„[y) ■ 



{y,t*;to) 



Jta JS.'^ Jta 

^^^xZ'..Ay), (B7) 
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which entails that ^ / {y,t;to) for any integrable f{y) is bounded in the time 

t dependence if and only if 

/ dyH°''^{y)f{y) = [ dy/(y) = . (B8) 
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